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We put forward a theory on the optical force exerted upon a dipolar particle by a stationary 
and ergodic partially coherent light field. We show through a rigorous analysis that the ensemble 
averaged electromagnetic force is given in terms of a partial gradient of the space variable diagonal 
elements of the coherence tensor. Further, by following this result we characterize the conservative 
and non-conservative components of this force. In addition, we establish the propagation law for 
the optical force in terms of the coherence function of light at a diffraction plane. This permits us 
to evaluate the effect of the degree of coherence on the force components by using the archetyp- 
ical configuration of Young's two apertures diffraction pattern, so often employed to characterize 
coherence of waves. 

PACS numbers: 



INTRODUCTION 

Since the works by Ashkin [l], Q , optical trapping and manipulation of particles became a tool of wide interest and 
use. Optical tweezers of micrometric objects are now frequently employed in various areas of science, particularly in 
biology. 

The mechanical action of light on particles is a consequence of its momentum [3|. In this context some studies 
concerning momentum conservation laws have been established both for deterministic 0-Q and partially coherent 
wavefields [1,0]. However, in spite of the vaste research on coherence of light (even in the context of inverse problems 
13- HI) and on optical trapping, no study on optical forces from fluctuating partially coherent fields has been developed 



yet, to our knowledge; apart from a couple of works -hj, [15( on Rayleigh particles in some particular focusing 
configurations. 



In the current knowledge on photonic trapping 16 21|, the illuminating light is assumed to be coherent, however, 
any field produced by a finite source is partially coherent [22T-|26j . Addressing the partial coherence of light in optical 
manipulation of objects should become increasingly important as one enters in the nanoscale (or subwavelength region) 
in near-field studies [l6| , and as a wider use is done of antenna-like sources with partial fluctuations as well as of the 
partial coherence induced in thermal sources. 

In this paper we analyze in detail the different contributions to the optical forces exerted on a small particle by 
random stationary and ergodic partially coherent external fields [22| . We emphasize the influence of the degree of 
coherence on these forces. To this end, we address a system that, since early studies, has been paradigmatic to observe 



the nature of light and matter waves, as well as to characterize the degree of coherence of wavefields [2^| . This is the 



Young interference pattern from two small apertures of an opaque screen. Concerning our study, we shall consider 



this configuration as discussed in the classical work by Thompson and Wolf 27] that relates the observed visibility of 
the interference fringes with the estimated degree of coherence of the light at those two apertures. 

Hence, we shall establish a theory for the mean optical force on a dipolar particle, (understood as that whose electric 
and/or magnetic polarizability is due to the corresponding first electric and/or magnetic Mie coefficient dHIl). tws 
includes the limiting case of Rayleigh particles. Then we will study the dependence of this force gradient, scattering 
and curl components on the cross-spectral density of the fluctuating stationary wavefield. We pay a special attention 
to the transition to a scalar theory [3l[ which characterizes most experiments in Fourier optics that do not measure 
or observe depolarization effects. In this way, we present in Section 2 a rigorous formulation for the ensemble average 
of the electromagnetic force; proving that this is given in terms of a partial gradient of the space variable diagonal 
elements of the wavefield coherence tensor, (this tensor being electric and/or magnetic, depending on whether the 
particle responds to the electromagnetic wave with an electric and/or magnetic polarizability). Since this coherence 
function obeys the Zernike propagation law [22| . we prove that the optical force on a small particle from a partially 
coherent wavefield depends on the field coherence tensor, or cross-spectral density tensor, on a plane from which this 
wave has propagated. 

Further, in Section 3 we consider the classical two-apertures configuration put forward in Thompson and Wolf 



's work [2jJ. We perform calculations of the diffracted field, obtaining interesting interference pattern distributions 
for the different Cartesian components of the optical force exerted on a dipolar particle situated in a plane in the 
Fraunhofer region with respect to that of diffraction. Thus, we shall characterize the different magnitudes of these 



2 



conservative and non-conservative force Cartesian components, both due to the diffraction process and to the degree 
of coherence of the light in the aperture plane, as well as steming from the particle polarizability. 



AVERAGED OPTICAL FORCE FROM A PARTIALLY COHERENT WAVEFIELD 



We shall consider fluctuating time stationary and ergodic fields 22, 3l|. For a single realization whose real elec- 
tric and magnetic vectors are E( r ) (r, t) and B^ r ^ (r,t), respectively, at a space point r and time t, the frequency 
decomposition is [22 1 



/oo 
-OO 

/OO 
-oo 



The corresponding complex analytic signals are [22t |31| : 

/OO 
E(r,oj)e- iuJt duj, 
-OO 
pOO 

B(r,i) = / B{r,oj)e- iult duj, 



where those Fourier integrals should be considered in the sense of distribution theory. In addition, we have 

E(r,w) = EW (r,w) ,u> > 
,w < ' 



(1) 
(2) 

(3) 
(4) 

(5) 



B(r,w) = BW (r,w) ,w > 
,w < 



(6) 



E(r,t) 
B(r,t) 



E' r » (r,i)+ffiW (r,i) 
BM (r,t)+<BW ( r ,t) 



(7) 
(8) 



The superscripts (r) and (i) denote the real and imaginary parts, respectively. For each Cartesian component of 
these electric and magnetic vectors, they form a Hilbert transform pair in the t variable 22, 3l| . 

We shall now calculate the ensemble average of the force exerted by the random field on a dipolar particle, (under- 
stood in the sense mentioned in Section 1), over its different realizations: 



(F(r, t)) = l (p« (r, t) ■ V) E« (r, t) + \ x B W (r, t) 



(9) 



p( r ) is the real part of the dipole moment induced by the fluctuating incident wave on the particle. If a e denotes the 
particle electric polarizability, one has that 



p(r,i) = a e E(r,i) 



(10) 



Denoting: dp/dt = p, and due to a well-known property of the derivative of Hilbert transforms [22j, the real and 
imaginary parts of each Cartesian component of p are Hilbert transforms of each other in t. 
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Let us evaluate the fist term of Eq. ©. Taking Eqs. ([3]), dU), (UJ, (JSJ) and (fT0|) into account, omitting the explicit 
r, t dependence in the forthcoming notation; this is (cf. [22|. |23|. l31j): 

(r) fl c(r) 



((p i + J »j)a i (£? i + ^)) 
i 



= Lim — - / dt 



T^oo2T 



Pj(r,ui)djEi(r,u) 2 )) duiiduj 2 



+ 
+ 



Pj*(r,uii)djEi (r,w 2 )> duiduj 2 



(11) 



Where i,j = 1,2,3 and Einstein's convention of omitting the sum symbol on the repeated index j has been 

used. Using the properties of Dirac delta distribution: 2-k8{lu) = $™ e~ lut dt and 5(uj) — 5(—u>), the previous equation 
reduces to 



Mi eiM 



2tt 



+ 
+ 



g^' E ^ (r, wi, W2) <5 (wi + W2) dujiduj 2 



gf ,B \r,LUi,u> 2 ) S (uji — u) 2 ) du>iduj 2 



'If ,£) (r, wi.Wa) <5 (wi - w 2 ) dwidw 2 



Sfp 3 ,£ ) (r,wi, w 2 ) 5 (wi + w 2 ) dujiduj 2 



where the cross-spectral function g^ U ' V \r, ojx,lj 2 ) is 



9i U ' V) (r,wi,u 2 ) = Lim^^; (U :i (r,u) 1 )d j Vi(r,u: 2 )) , 



(12) 



(13) 



Uj(r,u>) and Vi(r,w) (j, j = 1,2,3) being the spectra of two analytic signal Cartesian components. 

On performing the uj 2 integration in Eq. (IT21 and taking into account that due to Eqs. - @ , and to Eq. (fT31) . 
one has that 



9i (r,wi,-wi) = 0, 
only the second and third terms of Eq. (fT2"|) arc different from zero. Thus finally 



(14) 



1 



Lim ^ {pj^w^djEi (r,wi)) dwi 



4tt5R 

J — OO 

4ttK / ^'^(r,^,^)^!, 



(15) 



where 5ft denotes the real part. 

On the other hand, writing as eyfe the antisymmetric Levi - Civita tensor, (i, j, fc = 1, 2, 3); the second term of Eq. 
® is: 



ieyfc^B^) = l eijk ((p j+ p*)(B k + B* k )) 



(16) 
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C T->oo2T 



dt 



-T 



----- / !-^ 1 )e- l(wi ^ 2)t (^(r, Wl )B fe *(r, W2 ))rf Wl rf W2 
+ / iw 1 e+^ 1+W2 )*^--*(r,wi) J B fc *(r,a; 2 )) C L; 1 cL; 2 



(17) 



Or in a more compact form 



^jfc((Pi+Pi)(B*+Bit)) 
2^ r r°° 

(-iuj l )Wll ,B \r,uj 1 ,uj2)S(uj 1 + u; 2 )du>ida/2 



+ / (-iwi)W^ ,B '(r,wi,£i;2)5(£*;i -o^dwidwg 

/oo 
-00 

/oo 
mjiW^ ,B (r,cJi,o;2)^(^i +£^2)^1^; 
-00 

where the cross-frequency density tensor W^ (r ^1,(^2) is 



(18) 



W 



(19) 



On performing the w 2 integration in Eq. (fT8|) and taking into account that due to Eqs. (|4"| - (|6]l. one has that 



W 



(U,V) 



(r, wi, — wi) — 



(20) 



only the second and third terms of Eq. (|18j) remain different from zero. 

Now, since B = c/iuj'V x E, i.e. —iujBl — cekimdiE^, and taking into account that e%jh^hlm = o~no~jm — Si m Sji, Su 
being the Kronecker delta unit tensor, Eq. (|18l) becomes 



1 



c 

= 4tt3? 



Lim — — 

T-Kx>2T 



Pj(r,LJi)diEj (t,uii) 



Pj(T,U)i)djEi (T,U)i) 



duii. 



We introduce the electric field coherence tensor [22|,|23( £,fc(r,r',r) = (E 3 ;(r, t)E k *{v', t + r)) expressed as 



£ ife (r,r',T) = / £ jk (r,r',uj)t 



Sjk(r, r',u>) being the electric field cross- spectral density tensor defined as 

£ jk (r,r',u) = Lim — <^(r, u x )E k {T' , u 2 )) ■ 



(21) 



(22) 



(23) 
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On introducing Eqs. (fT5)l and (|2Tj) into Eq. ([9]), and taking Eqs. (fTUj) . (|22|) and (l23l) into acount, one finally obtains 
for the averaged force acting on a dipolar particle Q , the following expression in terms of the analytic signal associated 
to the random field 



(Fi(r,t)> - 4tt» / Lim^(^(r,^)9^/(r,^))da; 



— oo 
oo 



= 4kM I n e d't' ) Lim^ f {E s (v..:)E l (v.^VuU- 



= 47rK / aed^TrEjkir^w)^, (24) 



oo 



where we have substituted u>i by uj. The symbol Tr denotes the trace. On the other hand, d\*^ means that the 
derivative with respect to the ith component of r is made on the complex-conjugated component E* . 
Eq. (|24[) may also be expressed in terms of the coherence tensor as 

(F t (T,t)) = 47r5R{( Pj (r,i)a^;(r,t))} =4 7 rK{a e 9f ) <^(r,i)£;;(r,t))} 

= 4 7 rSR{a e 9|* ) Tr%(r,r,0)}. (25) 

In Eq. (|2l)|) we have recalled that acording to (|2"2"|) one has that £jk(r,r, 0) = f_ £jk{r,r,Lj)duj, On introducing the 
mean force Fourier components < Fi(r, w) ) as 



(Fi(r,t)> =2tt / ^(r,w))dw, (26) 

we see that according to Eq. (f^E)) we may express them as: 

'F<(r,w)\ = 2^[a e d ( l * ) Tr£ jk (r,r,uj)}. (27) 



The above calculations also lead to the conclusion that , being p, p, E and B analytic signals of t, one also has 
that 

p<f>B<p) = 2$t{p jB * k ), 



pV'diE?') = 2tf /-.y^/-:; . (28) 



Notice that introducing Eqs. (|2"B"j) into Eq. ©, one obtains again (|2"5| . 

Eq. (f2"5)) shows that the mean force is linked to the coherence tensor of the field. This latter quantity fulfills the 
Helmholtz equation whose integral representation leads to well known propagation laws in coherence theory, like the 
Zernike law and the Van Cittert-Zernike theorem, or to the dependence of the intensity on the degree of coherence of 
the wavefield in the primary or secondary source surface that emits it 22J, |23|, |3l| . Hence all these phenomena have 
consequences for the averaged force. 



Conservative and non-conservative components of the averaged optical force. The case of magnetodielectric 

particles 



It is well-known that the time averaged force from coherent fields may be expressed as the sum of three parts, one 
conservative and two non- conservative, (cf. 2^, 28[), namely, a gradient, a scattering and a curl of a electric spin 
density. Similarly, the force spectral components given by Eq. (|27[) lead to 



F(r,t)) = 27r5fta 

+- 47r3a e 5R 
f 47rQ ; Q;eQ ! 



+00 1 

Lim — ( V 

-oo T— >oo 2T 

+ OO 



E(r,w) 



2\ 



-00 
+00 



Lim — (fcE(r, u) x B* (r, u) ) doj 

T— >oo 11 \ 



T ^((E*(r >W )-v)E(r,a,) 



duj 



(29) 
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where 3 denotes the imaginary part. In Eq. (|29[) the first term represents the conservative or gradient force, whereas 
the second and third terms correspond to the non-conservative scattering component, or radiation pressure, and to the 
curl force, respectively. Likewise, one may write the same decomposition for the averaged force spectral components 
in bj— space, [cf. Eq. ([27)) ]: 



F(r,w)> = Ka e V( E(r,w) 



2fc3a e 5R{^E(r,cj) x B* (r,w)U 



23a e 3 | ( (e* (r, u) ■ v) E (r, w)\ } 



(30) 



It should be remarked that if the particle is magnetodielectric, namely, if additionally it has a magnetic polarizability 
c*m [HI , then in a similar way as for Eq. (|30l) one obtains for the averaged force on the particle due to the magnetic 

field FT (r, t) = 2K{a m 0j* ) TrB ifc (r,r,O)}, B jk (r,r',r) = (^-(r, t)B* k (v,t + t)>: 



F m (r, W )\ = »a m V 



B(r,w) 2 \ +2fcQa m 5ft{^E(r,w) x B* (r,w)J}} 



+ 



23a m 3 { ( (B* (r, w) • v) B (r, u)} } 



(31) 



And for the mean force due to the interaction between the electric and magnetic dipole induced in the particle 
(Fr m (r,t)) - -(8/3)k^{a e a* m e ljk g 3k (r,r,0)}, g jk (v,v',r) = (^(r, t)B* k (v, t + r)> 0: 



|sR(a e a^)5R^E x B*^> - S(a e <,)9^E x B*^| 



4fc 4 

IT 

- — 3t(a e a4)3?(E x B" 



4fc 3 



3(a e a4) 









[M 


E 











3t((E* ■ V)E 



(32) 



For the sake of brevity, we have omitted in the notation of Eq. (I3"2l the arguments r and <J of the analytic signal 
spectral vectors E and B. In this paper, we study the mean force on a particle with electric polarizability a e only, 
[cf. Eq. dSHJ)]- 



Dependence of the averaged optical force of propagated fields on the coherence at a diffraction plane. A 

Young interference configuration 



To illustrate the above with a simple example, let us consider a wavefield whose frequency components may be 
described by a scalar function U(y,uj). (This space-frequency description may also apply to a quasimonochromatic 
field, harmonically vibrating as exp(— iuit) around a mean frequency Q). In this case E(r, uS) — U(r,u)e(uj), (cf. 
Section 8.4 of [31]). This a common situation in Fourier optics [32]. The vector e(u) is real (linear polarization). On 

introducing W(ri,r 2 ,uj) — (U*(ri,u)U(r2,u))) as the cross-spectral density of U(r,u)) [22j], writing W(r, r, lu) = 
(U(r,uj)diU*(r,uj)) and taking real and imaginary parts in Eq. (|2"Tj) one obtains 

(Pifau)) = 2\e{uj)\ 2 ^{a e {d l U*{v,uj)U{r,uj))} 

= 2\e(uj)\ 2 [$la e $l(U(r,uj)d t U*(r,uj)) -$Sa e $i(U(r,u)diU*(r,Lj))]. (33) 
The first term of Eq. (|33p . is the mean gradient force, which is expressed as [see also the first term of Eq. (|3T))) ]: 



F? (r,w) 



= \e(u)\ 2 $ta e diW{r,r,w) 
= |eH| 2 3?a e a 4 (|C/(r, W )| 2 ). 



(34) 



The second term of Eq. ([33]) is proportional to the mean energy flow spectral density (S) associated to the scalar 
wavefunction U(r, w) [33j: 



(5i(r,w)> 



-Q{U(r,u)aiU*(T,u)). 



(35) 
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As such, it is the averaged scattering force, or mean radiation pressure, i.e. 

^ sc (r, W )) = -2|eH| 2 3a e 3{<9 4 ( *V(r,r, W )} 

= -2|e( W )| 2 3a e 3{([/(r,^)a 4 (7*(r, W ))} 
= 2k\e{oj)\ 2 Sa e (S l {r,L J )), 



(36) 



Eq. (|36p manifests the correspondence of (<S) in this scalar formulation of the radiation pressure with the mean 
Poynting vector (S) = (c/8tt) (E x B*) in the second term of Eq. (|3T)|) acording to the vector representation. 

Notice that the mean curl of electric spin density which according to the third term of Eq. (|30[) takes on the form 



F' 



curl 



(r,w)) - 2<3a e $s{e*(u>)e i (Lj)d j W(r,r,u>)} , 



(37) 



where di means the derivative in the non-conjugated wavefunction, i.e., diW(r,r,ui) = (U*(r,uj)diU(r,uj)). Now, 
taking into account Maxwell's divergence equation V • E(r, u>) = in Eq. (f3"T|). it is easy to demonstrate that 

F[ wrl (r,u>)^ will be zero if 9 {e*(u>)ei(u>)} = 0; which evidently holds since e(u>) is real. 

According to the Huygens-Fresnel principle 22, 23, 3l| the fluctuating field propagated from points r' of a surface 
A up to a point r is given by 

U(r,u,) = -^- J U(T>,u)^<Pr>. (38) 

where R = \r — r'| and k = 2tt/X, the wavelength being A . Thus from Eq. (|25|) the averaged force on a dipolar 
particle in r will be 



Fi(r,u) 



2[-) |e(u;)| 2 ^a c 



Ri R2 



[ ik+ ±- ] R ^W(r' 1 ,4,a;) 

AJA u l u l 



d 2 r[d 2 r' 2 } , 



(39) 



Ri = r — r'i, Ri = |r — rj|, (i = 1, 2). As mentioned above, Eq. (|39| exhibits the dependence of the mean force exerted 
by the propagated field on its coherence properties on a surface A. 

For instance, we consider the surface A being composed of an opaque screen with two point holes, (see Fig. 1), so 
that the random field wavefunction in A is: U(r',u!) = U(qi,ui)5(r' — qi) + f/(q2, u))S(r' — q%). Then from Eqs. (|38l) . 
(l39l) . (j34l and (j36j) we obtain for the conservative and non-conservative force components on a particle at a point P 
of position vector r: 



fgrad ^ 
2 



2tt 



|eH| 2 »a e { <|C/(qi,^)| 2 ) ^ + <|t/(q 2 ^)| 2 ) * 



R 2 



|W r (qi,q 2 ,w)| 


YRi 




R1R2 


.Uf 





cos (fc(Ri - R 2 ) + a(qi, q2,w)) 



R R 

— ^ - ) fcsin(fc(Ri - R 2 ) + a(qi,q 2 ,w)) 



(40) 



+ 



+ 



F sc (r J a;) / 
|W(qi,q 2 ,w)| 



Ri 



(kRi cos (fc(Ri - R 2 ) + a(qi, q 2 , w)) 



i?i-R 2 
fc(Ri - 

(fci? 2 cos (fc(Ri - R 2 ) + a(qi , q 2 , u>)) 



sin(fc(Ri - R 2 ) + a(qi,q 2 ,w))) 

R_2 

Rl 

sin(fc(Ri - R 2 ) + a(qi,q 2 , w)))]} 



(41) 
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FIG. 1: Schematics of the configuration for observing interference at points P = (x, y, d) of a screen B by diffraction of light, 
propagated from an incoherent source, at two apertures in A, centered at points r[ = (qi,0) and r' 2 — (q2,0), respectively, 
qi = {xi,yi), q2 = (3:2,2/2). 2h = |qi - q 2 |. 



Denoting Rj = r — qi, (i = 1, 2), and a (qi, q 2 , w) being the phase of W(qi, q 2 , w). 

In the Fresnel and Fraunhofer regions one may approximate R\ ~ i? 2 in the denominators of Eq. (|40j) . Also, for 
kRi 3> 1, (i = 1,2), the cos terms are negligible versus the sin terms and the gradient force has an interferencial 
sin behavior, proportional to the difference: Ri — R 2 . On the other hand, the sin terms of (|4ip are negligible 
versus the cos terms, rendering a scattering force proportional to the intensity pattern. This will be discussed again 
in Section 3 in connection with the configuration of Thompson and Wolf experiment, which replaces the two point 
holes of this schematic example by real apertures. In addition, by dropping in Eqs. (140)) and (|4ip the corresponding 
factor constituted by the real and imaginary part of the electric polarizability, we observe that the action on particles 
situated at points Ri 3> A by the SJa e -normalized repulsive scattering force produced by each independent pinhole, is 
much larger along R; than that of the corresponding 5fta e -normalized attractive gradient force. 



INTERFERENCE OF TWO RANDOM WAVES: DEGREE OF COHERENCE AND AVERAGED 

OPTICAL FORCE 

In this section we address the force on a dipolar particle in the configuration of the classical two-aperture arrange- 
ment by Thompson and Wolf, employed in 1957 to observe and characterize the degree of coherence of a light wave 



27]. As shown in the scheme of Fig. 1, the wavefield emitted by an incoherent source is brought by a lens L\ to the 
mask A in z = 0, containing two small circular apertures of radius a, centered at points r[ = (qi, 0) and r' 2 = (q 2 , 0), 
respectively. A second lens L 2 sends the field diffracted in A to points P of a screen B coinciding with its focal plane. 
The degree of coherence of the wave in A: MisU i92> u ) = W(qi, q 2 , w)/[W(qi, qi, w)PF(q 2 , q 2 , w)] 1 / 2 is expressed by 
means of the Van-Cittert-Cernike theorem [22j, [31| in terms of the intensity exiting the incoherent source . In this 
way, we stablish the influence of the partial coherence in A of the light emitted by the random incoherent source, on 
the optical force from the diffracted field upon a dipolar particle placed in B. 

In the vector theory of diffraction, within the range of validity of the Kirchhoff approximation, the diffracted electric 
vector produced by an aperture centered in r = in a screen A is expressed in the far zone as 0, [33[ 



E(r,w) = ^— kx f nxEW(r»e- ,kT 'ds', (42) 
27rr- J A 

where k = fcs = (27r/A)s, s = r/r is a unit vector in the direction of observation r = (x, y, z), r' denotes a coordinate 
in the aperture whose element of surface area is ds' , and n is the unit outward normal to ds' . The time-dependence 
e -iwt j s understood, and E^^(r') = (w) exp(ifcnW • r'), (|n^| = 1, 3?eW = 0), is the electric field incident on the 
mask A. 
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Assuming n.W = (0,0,1), a flat opaque screen in z = with a circular aperture of center r' 
produces according to Eq. (|42|) the diffracted field [3l], 



E(r,w) = U(r,u)e(u)) 



and radius a 



(43) 



where 



E7(r,w) 



7ra 



f'o 



(44) 



In Eqs. (|43|) and (|44l) we have written e(w) = s x (n x e^-'(w)) and vq = ka-^x 2 + y 2 /z, respectively. 

Next we study the effect of the degree of coherence in A of the waveheld propagated from the chaotic source on the 
mean force upon a dipolar particle in the Fraunhofer zone. Hence we shall address the consequences of the correlation 
between the held at the two circular apertures in z = 0, centered at qi = (0,h,0) and q2 = (0, — h,0), as shown in 
Fig. 1. To this end, we evaluate the field at an arbitrary point P in the plane B at z = d where the particle is 
situated, produced on diffraction in A, (cf. Fig. 1) within the Kirchhoff approximation: 



U(v,oj) 



ikd 

e e 



iXd 



U(qi,u)e 



U(q 2 ,uj)e i 



(45) 



where U{(\i, ui) (i = 1, 2) is the complex amplitude of the random wavefield at qi and q2 emitted by the fluctuating 



source. If as in Thompson and Wolf experiment [27], one has that ( |f/(qi,u;)| 
mean intensity at P then is known to be 



, Ji Oo) 



In, the observed 



(I(r,u>)) = (E(r,w) -E*(r,w) 



= 2/ 



ira 2 |e (uj)\ 
Xd 



11/ \ i ( , / \ 2khy\ 
1 + |/4qi,q 2 ,a;)| cos I 0(q l5 q 2 ,w) H ^— I 



(46) 



The factor 2khy/d represents the path difference |Ri — R 2 |, (see Fig. 1). 0(qi,q 2 ,w) is the phase of the degree 
of coherence: /x(qx,q 2) w) = 2J\{u)/u in z — 0. Where u = 2nph/(XA). p being the radius of the source, assumed 
planar and circular, and A denoting the distance between the source and L\ 27 1. 



The interference law of Eq. (|4^]) is well known [27], |3l[ . We shall perform calculations of the force with the same 
parameters as in Ref. [27j ]. namely: A = 579 nm, 2h — 6mm, a — 0.7mm, d = 1.5m. We consider 0(qi,q 2 ,w) = 
and we will initially normalize the results to the polarizability, so that we will make 3?a e = 3a e = 1; this allows us 
to obtain an estimation of the relative strengths of the different force components due to diffraction, independently 
of the polarizability. 

In this far zone, the gradient force is governed by the expression (|34[) applied to the mean intensity (1461) . Since 
the apertures are aligned along the y— axis, the y component for kRi — > oo (i — 1,2) is obtained after a long but 
straightforward algebra 



pt-i 
y 



rad 



-45fta e /o 



wa 2 \e(ui)\ 



Xd 

x |/i(qi,q 2 ,w)|sin 



hk 
~d 



. 2khy\ 
(qi,q 2 ,w) + — — I 



(47) 



which agrees with the remark at the end of Section 2 concerning Eq. (|40)) . 

This expression is just the derivative of Eq. (|46|) with respect to y, assuming that the factor outside the brackets 
in (|46| is constant, (although this is not strictly true, the terms yielded by the y-derivative of this factor become 
negligible, as shown in the Appendix A). The other two components: (^F^ rad ^ and (^F9 rad ' K j arG similarly obtained in 
the Appendix A. 

Fig. 2 shows the interference pattern of (/), normalized to its maximum, at the screen B for |/x(qi,q 2 ,w)| = 1, as 
well as the spatial distribution of the three components of the mean gradient force on a dipolar particle in B due to 
this distribution of light. To see the relative weight of each Cartesian component, we normalize it to the magnitude 
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FIG. 2: Spatial distributions in the XY plane of the intensity and the normalized gradient force components for |/x(qi, q2, w)| = 
1. (a) Normalized mean intensity (J), (b) ^F| rod Y (c) (Pjj rad \. (d) {Ff 1 ^. All values are calculated on a dipolar particle 
at the screen plane B, placed at distance z = d — 1.5 m from the aperture screen A. The force components are normalized to 
Ka e and to the magnitude of the total mean force |(F' ot \| = ]/F 9rQd \ + (f sc \J. 



of the total mean force 



■pgrad 



We also observe an interference pattern along OY in each 
component of this conservative force, (^F^ rad ^ being much smaller than the other two. In addition, Fig. 2(b) exhibits 
an oscillatory modulation of \F^ rad \ along OX, (cf. Appendix A). 



We remark that in the limiting case |/x(qi, q2, w)| = 0, {F~j rad j is just proportional to d y (2J\ (vq) /vq) 2 and the 
interference effect disappears, as it should. Since this y derivative was neglected versus the term kept in Eq. (|47p . 
(see also Appendix A), the values of I^F^ rad ^j then are practically zero compared to those due to a partially coherent 
wave. This is seen in Fig. 3. The intensity pattern, which acts as a potential distribution for the illuminated particle, 
is shifted by 7r/2 with respect to that of the conservative force (^F~j rad ^j, whose oscillation amplitude progressively 

diminishes to zero as the value of |/x(qi, q2, w)| decreases. This behavior of the conservative force constitutes the 
basic mechanism of an optical tweezer with several equilibrium positions of the particle along the lines in the screen 
B where (I) is maximum, . Such points occur along OY at x = , [cf. Figs. 2(a) and 2(c)], and are precisely those 

where ^F| rad ^> = 0, [cf. Fig. 2(b)]. 

We next address the scattering force on a small particle in an arbitrary point of the screen B, obtained on introducing 
Eq. (|45p into Eq. (|36l) . (see Appendix B). The scattering and gradient force x— components are of similar magnitude, 
but of signs opposite to each other; this is seen on comparing Fig. 4 (a) with Fig. 2 (b). By contrast, the y-component 

of ^F sc ^j is one order of magnitude smaller than its homologous of the gradient force, suffering a sharp change of sign 
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y (mm) 

FIG. 3: (a) Normalized mean intensity (I), (b) Normalized gradient force component (^Fy rad ^ for different values of 
|/x(qi,q 2 ,w)|. 



at y — 0. However, ( F 2 SC > which is given by (see the Appendix B) 



4fcw ( W|eH| y(V iM 



Arf 



/ 2khy\ 
1 + |/i(qi,q 2 ,w)| cos I 0(qi,q2,w) H — I 

2fcSa e (/(r,w)), 



(48) 



is seven orders of magnitude larger than the corresponding conservative force i^F^ rad j, [compare Fig. 4 (c) with Fig. 

2 (d)]. This, which is in accordance with the remark of the last paragraph of Section 2 concerning Eqs. (l40l) and (|4"Tj) 
for waves from two pinholes, stems from the proportionality of (F£ c ) to the Poynting vector [28[ and hence to the 
mean scattered intensity in the far zone; [observe that the normalized (F^ c ) of Fig. 4(c) is identical to the normalized 
mean intensity (/} of Fig. 2 (a)]. As a consequence of the conservation of momentum, the particle is pushed towards 
z > d (see Fig. 1). The ratio between the maximum values of the gradient and scattering force components, (see Eqs. 

dH]), (08J and the Appendix A and B), /Pf ad ^ / = -(y/d) (P<> rad \ / (f° c ^ = (yh/d 2 ) /l + |/z|) explains 

the difference between the magnitudes of these force components. 

The much larger strength of the normalized scattering force may prevent the lateral manipulation of the particle 
in B. If this were the case, it may be overcame with a scheme analogous to that employed in holographic optical 
tweezers 



34 



Summing up, we observe that the Young experiment configuration shows us fundamental characteristics of the 
optical force components, which in this system allow a scalar formulation and thus yield no curl component. The 
mean scattering force is proportional to the mean scattered intensity, its longitudinal z-component being several orders 
of magnitude larger than its x and y components. On the other hand, the y and z-components of the gradient force 
are proportional to the magnitude of the degree of coherence, thus becoming zero for incoherent light, and hence their 
value oscillates and decreases as Ji(27rp/i/AA)/ (2irph/ XA) versus the distance h between apertures [27j | . 



Effect of the electric polarizability on the mean optical force 



We have so far estimated the different components of the mean force by only considering the configuration of the 
diffracted waves; namely, by normalizing them to the particle polarizability. However, it is worth remarking that since 
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FIG. 4: Spatial distribution in the XY plane of the normalized scattering force components for |/x(qi, q2, = 1. The 
normalization factor is 3a e ^F*°*\ . (a) (F^ c ). (b) (Fy C ). (c) (F^ c ). All values are calculated on a dipolar particle in the 
screen plane B, placed at distance z = d — 1.5 m from the aperture screen A. 



their actual strengths, observed in an experiment, are proportional to ^fta e (gradient force) ant to Ssa e (scattering 
force), the relative values of these two parts of a e should greatly influence the magnitude of these forces. Notice 



that although we have concluded that (F^ r j <C \ Fz°j when they are normalized to J?a e and 3a e , respectively, in 

most cases pertaining to dielectric particles one has that SRa e 3> 5a e , except in the presence of Mie electric and/or 
magnetic 28, 3^] or plasmon j35| resonances. 

For a small spherical particle of radius r*o, with relative permittivity e p , in the Rayleigh limit (kro <C 1), we adopt 
the expression for the dynamic electric polarizability [28j, [30l ] : 



i-k 3 a^ 



being the static polarizabilty 



„3 C P 



(49) 



(50) 



As an illustration, we consider a dielectric particle with ro = 25 nm and e p = 2.25. With these data, we observe 
as mentioned above that $ta e = 4593 nm 3 ;§> 3a e = 17 nm 3 . The illumination that reaches each aperture is assumed 
with a magnitude of the Poynting vector (c/2)Iq \e(uj)\ 2 — 10 12 W/m 2 . Fig. 5 shows the corresponding different 
components of the total force, (this time of course without performing any normalization). 

As seen, the patterns of Fig. 5 (a) and Fig. 5 (b) are equal to those of Fig. 2 (b) and Fig. 2 (c) respectively; this 
implying that the scattering force is negligible compared to the gradient force along OX and OY. However, although 
5fta e ^> Sa e , this is not enough for (^F^ rad ^ to exceed (f^°^, (remember that we obtained a difference of seven 

orders of magnitude between these two normalized z-components), therefore the contribution of /fijj ra <^ to (P^ ot ^ i s 
negligible by four orders of magnitude. Notwithstanding, it is important for trapping purposes that the y component 
of the force (^Fy 0t \, which is of conservative nature, is of the same order of magnitude as the non-conservative z-fovce 

F* ot ^. As the coherence diminishes, Fig. 3(b) gives an assessment of the corresponding decrease to be expected in 



both 



' 'F* ot ^ and (^Fy 0t ^ from their values in Figs. 5(a) and 5(b). 
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FIG. 5: Spatial distribution, in pN, of the total force Cartesian components on a dielectric particle with ro = 25 nm and 
e p = 2.25, in the screen plane B, placed at distance z — d = 1.5m from the aperture mask A, |/x(qi, q2, w)| = 1. No 
normalization is done, (a) (F* ot ). (b) (F* ot ). (c) (Fj ot ). 



CONCLUSION 



We have established a theory of the averaged optical force exerted by a partially coherent stationary and ergodic 
electromagnetic field on a dipolar particle. Although we have put forward expressions for the force acting on magne- 
todiclectric particles, we have illustrated the effect of partial coherence on a particle that has electric polarizability 
only. This has been carried out by calculating the influence that either the cross-spectral density and the degree 
of coherence of the fluctuating field at a certain plane of propagation, have on the force components due to light 
propagated beyond that plane. 

To this end, we have made use of the classical two-apertures interference configuration of Young to calculate the 
force on a dipolar particle in the Fraunhofer zone. Hence, in correspondence with the classical study on coherence by 
Thompson and Wolf [27| , we have calculated the force components and their dependence on the degree of coherence 
of the fluctuating field at the plane of the apertures. The result is quite interesting; it shows a fringe pattern spatial 
distribution in each Cartesian component of the conservative and non-conservative forces, and illustrates the way in 
which they compete with each other to the resulting total forces. 

The case of magnetodielectric particles may be similarly studied. Then, apart from adding to the above the pure 
magnetic force whose strength depends on the magnetic polarizability of the particle, an electric-magnetic interference 
force component, of opposite sign to the former, has to be summed to the latter. The relative weight of these latter two 
forces, as well as their dependence on the degree of coherence of the fluctuating field, should be evaluated. Of special 
interest will be the case of nanometric size high index, or semiconductor, spheres, illuminated in the near-infrared at 
which both their electric and magnetic dipole Mie resonances are excited [29|, |30( . 

Likewise, although we have emphasized our illustration of Section 3 with an scalar theory, it will be of interest to 
calculate these optical forces in terms of the degree of polarization of the fluctuating fields by making use of the full 
electromagnetic model of Section 2. 

We believe that this study will motivate further research on particle manipulation by fluctuating wavefields, both in 
the context of partially coherent waves and in the more general of random wavefields 36|- 3t| , which should permit to 
work with light and other electromagnetic wave sources like partially fluctuating nanoantennas of limited coherence 
or correlation length or thermal sources in which partial coherence is undeccd. This will be of special importance at 
the nanometric scale in the near field, where new consequences are to be found. 



Appendix 



Analytical expressions of the gradient force 
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rad 
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v = ka^/x 2 + y 2 /z. Observe that in this equation, none of the two terms may be neglected. 
Calculation of Ff ad 
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Calculation of Ff ad 



Therefore 



Calculation of F Q 



Calculation of Fy 



Calculation of 
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Analytical expressions of the scattering force 
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